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CALCULATION  OF  ATMOSPHERIC  COOLING  RATES  USING  MODTRAN2 


Lucia  M.  Kimball 
Graduaca  Student 
Mathematics  Department 
Worcester  Polytechnic  Institute 


Abstract 

A  technique  tor  calculating  atmospheric-  cooling  rates  has  been 
developed  based  on  infrared  radiance  calculations  using  MODTRAN2 . 
Comparisons  with  benchmark  line-by-line  calculations  show  very  good 
agreement.  The  technique  provides  significant  computational  time 
savings  associated  with  band  models  vs.  line-by-line  calculations. 
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CALCULATION  OF  ATMOSPHERIC  COOLING  RATES  USING  MODTRAN2 


Lucia  M.  Kimball 


Introduction 

Accurate  calculations  of  the  transfer  of  radiation  through  the 
atmosphere  are  necessary  for  many  climate  programs.  Line-by-line 
models  such  as  F.ASCODE  perform  detailed  computations  and  are 
considered  "exact"  within  the  limitations  of  the  current  knowledge 
of  molecular  properties  such  as  lineshape,  line  mixing  and 
interaction  with  foreign  gases.  However,  the  amount  of  time 
required  for  most  calculations  is  prohibitive.  M0DTRAN2,  the  most 
current  version  of  MODTRA.N,  the  Moderate  Resolution  Atmospheric 
Radiance  and  Transmittance  Model,  is  a  band  model  based  code  which 
has  been  extensively  validated  against  FASCODE.  For  a  500  cm'* 
spectral  interval  with  comparable  vertical  layering  the  time 
improvement  is  more  than  a  factor  of  100. 

In  this  paper,  we  consider  the  use  of  M0DTRAN2  to  calculate 
clear  sky,  infrared  cooling  rates  for  several  model  atmospheres. 
The  cooling  rate  results  are  compared  to  similar  calculations 
performed  with  FASCODE  as  developed  by  Clough,  st  al .  ,  (1992)  .  The 
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current  method  follows  that  of  Clough  as  closely  as  possible  so 
that  differences  in  the  results  can  be  attributed  primarily  to  the 
differences  between  the  line-by-line  and  band  models. 

Radiance  Calculation 


The  radiant  energy  emitted  at  wavenumber  v  along  a  path  with 
optical  depth  and  transmittance  is  given  by 

I,=]\'3jd(T!))dT', 

I^  =  f  ’3^  (d  (%)  )  . 

J  0 


where  3^  is  the  Planck  function  at  temperature  9, 


a, (9)  = 


2h 

gv 

e  -1 


Here  h  is  Planck's  constant,  c  is  the  speed  of  light  and  k  is 
Boltzmann's  constant.  In  the  present  application,  M0DTRAN2,  a 
moderate  resolution  transmittance  code,  has  been  used  to  first 
calculate  the  path  transmittances  and  secondarily,  the  radiances. 

The  version  of  MODTRAN  used  to  perform  the  necessary 
computations  utilizes  a  one  term  Pade  approximate  for  the 
effective  Planck  function  of  the  form 

3^  (aTj  5^  (9jj)  ]  (l^at  J 
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where  9^  is  the  temperature  ae  the  nearest:  boundary  of  the  Layer 

(the  upper  boundary  for  downwelling  radiance  and  the  Lower  boundary 
for  upweLLing  radiance)  ,  3,,=3(H)  where  H  is  the  mean  temperature 
for  the  Layer  1.  The  vaLue  chosen  for  a  is  0.2  after  Clough,  et 

al . ,  (1992).  For  a  band  modei,  is  an  effective  opticaL  depth 
derived  from  the  Log  ratios  of  the  fuLL  path  transmittances  between 
adjacent  Layers 


t.,  =  Ld 


The  radiance  for  a  singLe  Layer  is  then  given  by 


r  =  ( T’  -'T  )  3  (-  ) 

~v;l  '  -v.J  -v.l-l'  '  -v'  • 

The  fuLL  path  radiance  for  a  path  encompassing  muLtipLe  Layers 
is  found  recursiveLy  where  the  vaLue  at  each  Layer  depends  on  the 
contribution  from  the  previous  Layer,  see  Clough,  et  al . ,  (1981). 
Treating  the  ground  as  a  biack  body,  the  boundary  condition  is 
given  by 


/ 

V 


0 


Comparison  of  Radiance  Calculations  FASCODB  vs.  MODTRAN2 

The  most  accurate  technique  -  for  calculating  atmospheric 
transmittance  and  radiances  employs  a  Line-by-Line  code  such  as 
F.ASCODE.  Calculations  are  performed  at  small  enough  wavenumber 
intervals  to  adequately  represent  the  location,  strength  and  shape 
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of  each  spectral  Line  within  a  layer,  thus  justifying  Beer's  law 
for  calculating  the  optical  depth.  A  more  detailed  description  is 
available,  in  ELlincs::'.,  et  al.,  (1991).  Due  to  the  large  n:-.-’-"--  ■■ 
spectral  lines  this  is  a  slow  and  cumbersome  process  which  is 
therefore  impractical  for  most  uses. 

M.0DTRAN2  employs  a  two  parameter  band  model  to  offer  more 
efficient  computations  of  the  line  contributions.  Treatment  of 
line  tails  and  continue  is  similar  for  both  models.  The  two 
parameters  are  an  absorption  coefficient  (S/d)  and  a  line  density 
parameter  (1/d)  given  by 

(5/d)  =  (-^) 

Av  ^ 

where  Ay  =  1  cm"'  and 


(1/d)  =  (2./S7)  VSS... 

Av  *  ■ 

S;  is  the  integrated  line  strength  of  line  i,  and  the  sum  is  over 
N  lines  whose  centers  are  contained  in  a  bin  of  width  Ay.  The 
transmittance,  T,  for  the  spectral  bin  Ay  is  given  by 


,  -5  u  i>(v) 


Here  b(y)  is  the  Voigt  line  shape  function,  u  is  the  absorber 
amount,  S  and  n  are  defined  by  the  band  model  parameters  as 
follows: 


(s/d) 

(1/d) 
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n=  (1/ d)  4v  . 


For  a  more  detailed  description  see  A.  Berk,  et  al . ,  (1987). 

In  one  specific  comparison,  see  L.  Aireu,  et  al . ,  (1993), 
M0DTRAM2  was  shown  to  be  100  times  faster  than  a  non-optimized 
version  of  F.^SCODE  for  a  500  cm"'  spectral  interval  with  comparable 
vertical  layering.  With  proper  vector ization  and  parallel 
processing,  FASCODE's  speed  can  be  greatly  enhanced.  However, 
MODTRAN's  structure  is  equally  vectorizable ,  so  the  ratio  of 
improvement  should  ultimately  be  maintained. 

Flux  Calculation 

The  upwelling  and  downwelling  fluxes  of  radiant  energy  at  a 
given  atmospheric  level  are  calculated  by  integrating  the 

radiance,  ,  at  wavenumber  u  over  the  appropriate  hemisphere, 

Jo  JO 

where  jn  is  the  zenith  direction  cosine  and  is  the  azimuthal 
angle.  In  our  case,  it  is  assumed  the  radiances  are 
azimuthally  independent;  therefore  the  integrals  become 

(p)pG!JL  . 

J  0 
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The  required  integracion  for  determining  the  upwelling  and 
downwelling  fluxes  was  computed  using  a  standard  two  point,  first 
moment  Gaussian  quadrature.  Comparison  of  two  point  vs.  three 
pair.:,  ruadratur  a  ■  i "  ;t  show  sufficient  improvement  to  justify  tha 
additional  time  involved  in  the  computation.  When  simultaneously 
plotted,  the  2  angle  and  3  angle  cooling  rate  curves  were  nearly 
coincident.  The  direction  cosine  values,  associated  angles  and 
quadrature  weights  are  given  in  the  following  table: 


Direction  Cosine 

Angle 

i 

Weight  i 

0 .3550510257 

69.2034233 

0.1319536133  | 

i 

1 

0 . 3449489743 

32.3335307 

0 .3130413817 

Table  1:  Gaussian  quadrature  angles  and  weights 


MODTRAN2  was  used  to  calculate  the  radiance  along  a  path  to  space 
corresponding  to  each  of  the  given  angles  at  each  defined 
atmospheric  level.  The  downwelling  flux  was  then  found  by  the  first 
moment  quadrature.  The  upwelling  flux  resulted  from  integrating 
the  radiances  along  paths  corresponding  to  the  complement  of  each 
given  angle. 

The  net  flux  at  any  level  was  then  calculated  as  the  difference 
between  the  upwelling  and  downwelling  fluxes. 


Cooling  Rate  Calculation 

The  divergence  of  the  net  flux  at  atmospheric  level  I  represents 
the  rata  of  energy  los?  oo"  unit  volv- ^  --  atmosphere,  or  the 
cooling  rate 

Q  =dF  . 

The  monochromatic  cooling  rata  for  the  atmospheric  layer  bounded  by 
levels  I  and  l-l  is  computed  in  terms  of  change  in  temperature,  0, 
with  respect  to  time,  t,  by  the  finite  difference  formula 

00  I  _ 

where  g  is  gravitational  acceleration,  is  the  specific  heat  of 

air  at  constant  pressure,  and  P  is  pressure.  A  value  of  B.422 

g 

(mbar  K  d"')  per  (W  m'*)  is  used  for  the  ratio  ~  as  in  Ridgvay  et 
al . ,  (1991),  independent  of  altitude. 

Comparison  of  Cooling  Rate  Results 

All  M0DTRAN2  calculations  were  based  on  60  atmospheric  layers  as 
defined  by  Clough,  et  al . ,  (1992).  Layers  were  spaced  at 

increments  of  approximately  20  mb  pressure  from  0  mb  to  1013  mb. 
Cooling  rates  were  found  to  be  sensitive  to  the  chosen  layering  at 
both  the  top  and  bottom  of  the  defined  atmosphere. 
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Data  for  layer  temperatures  and  water  vapor  densities  were 
interpolated  from  the  published  ICRCCM  data  in  Ellingson,  et  al . , 
(1991)  . 

The  cooling  rate  results  as  compared  to  the  published  FASCODE 
results  were  very  good.  .Agreement  was  within  the  2%  documented 
FASCODE  -  MODTRAM  agreement  in  most  regions.  Relative  errors  in 
the  net  fluxes,  mid-latitude  summer  atmosphere,  were  all  less  than 
1%  with  the  exception  of  the  net  flux  in  the  ground  layer  which  had 
the  maximum  relative  error  of  3.7%. 

comparison  of  cooling  rate  results  for  mid-latitude  summer, 
tropical,  and  sub-arctic  winter  atmospheres  is  given  in  figures  1- 
3.  Plots  represent  spectrally  integrated  cooling  rates. 

A  representation  of  spectrally  dependent  cooling  rates 
corresponding  to  published  FASCODE  results  {Clough,  et.  ai . , 
(1992))  is  given  in  figure  4. 
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COOLING  RATE 


Comparisca  of  scecoraLly  Lnoegrated  cooling  races: 
ICRCCM  nid-'_acicude  summer  acmosphere,  water  vapor 
only.  Solid  line  represents  H0DTRAN2  results  and 
dashed  line  represents  FASCOOE  results. 


COOLING  RATE 

Comparison  of  spectrally  integrated  cooling  rates 
ICRCCM  sub-arctic  winter  atmosphere,  water  vapor 
only.  Solid  line  represents  MobTRAN2  results  and 
dashed  line  reoresents  F.ASCCDE  results. 


Figure  3:  Comparison  of  specPralLy  incegrated  cooling  rates: 

ICRCCM  tropical  atmosphere,  water  vapor  only.  Solid 
line  represents  MOOTRAN2  results  and  dashed  line 
represents  FA3C0DE  results. 
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